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Consider  the  boundary  value  problem  ty"  =  (y  -  t  )y‘ ,  -1  <  t  <  0, 

y(-1)  c  A,  y(0)  =  B.  Depending  on  the  choice  of  A  and  B,  one  can  insure 

*  *  2  "2 

the  existence  of  "turning  points"*  t;  y(t,£)  -  t  *=  0.  However,  due  to  the 
nonlinear  nature  of  the  problem,  one  does  not  know  the  position  or  number  of 
such  turning  points.  In  the  case  when  A  >  0  -  B  Kedem,  Parter  and 
Steuerwalt  gave  a  development  of  this  problem  based  on  an  abstract  bifurcation 
analysis  which  in  turn  was  based  on  ‘‘degree  theory^.  In  this  paper  we  give  a 
complete  analysis  of  the  problem  based  entirely  on  a-priori  estimates  and  the 
^shooting*^  method. 


AMS  (MOS)  Subject  Classifications;  34B15,  34E15 

Key  Words:  Ordinary  differential  equations.  Singular  perturbations.  Turning 
points.  Boundary  layers,  Multiplicity  of  solutions.  Asymptotic 
behavior 

Work  Unit  Number  1  (Applied  Analysis) 


*Will  also  appear  as  University  of  Wisconsin-Madison,  Department 
of  Mathematics,  Technical  Summary  Report  s2165. 

^Department  of  Mathematics,  Northeastern  University,  Boston, 

MA  02115. 

(2) 

Department  or  Mathematics  and  Computer  Sciences  Department, 
University  of  Wisconsin-Madison,  Wisconsin,  WI  53706. 


Sponsored  by  the  United  States  Army  under  Contract 

No.  QA,u329  -dO-C-dOl  L ,  by  the  Office  of  Naval  Research  under  Contract 
No.  N0001'i-76-0-0341,  and  by  the  National  Science  Foundation  under 
Grant  NSF-MCS77-0J715-01. 


A  COMPLETE  ANALYSIS  OK  A  MODEL  NONLINEAR  SINGULAR  PERTURLAT ION  PROBLEM 


HAVING  A  CONTINUOUS  LOCUS  OF  SINGULAR  POINTS  * 

Nancy  Kopcll^  '  and  Seymour  V.  Parter*^ 

1.  Introduction 

The  nonlinear  boundary  value  problem 

Ey"(t,e)  =  ty2(t,e)  -  t2]y’(t,e)  ,  (1.1) 

y(-1,£)  =  a,  y(0,c)  =  B  ,  (1.2) 

6  >  0,  was  introduced  by  Howes  and  Barter  [21  as  a  model  problem  having  a  continuous 

*  2  *  “2 

locus  of  potential  "turning  points",  i.e.  points  t  at  which  y  (t,E)  -  t  =0.  The  main 
questions  raised  by  (1.1),  (1.2)  qpneern  the  multiplicity  of  solutions  for  a  given  A 
and  B,  and  the  asymptotic  behavior  of  these  solutions  as  E  ♦  0. 

In  [2) ,  Howes  and  Parter  showed  that,  if  0<B<A<1,  the  only  possible  constant 
limiting  solutions  are  .  y  =  A,  y  =  B  and  y  5  1//T.  A  further  analysis  was  carried  out  by 
Kedem,  Steuervalt  and  Parter  (3J,  who  studied  the  case  A  >  B  =  0.  This  work  was  motivated 
by  a  conjecture  of  Sutton  [5]  ,  who  has  made  computational  experiments  based  on  the  methods 
of  [2].  Her  conjecture  was:  Suppose  0  <  B  <  1//3  <  A  <  1.  Then  for  c  sufficiently 
small,  there  exist  at  least  three  solutions,  having  A,  B  and  1//3  respectively  as 
limiting  values  for  t  e  (-1,0).  (Howes  [1J  has  since  given  a  proof  of  the  conjecture.) 

In  [3]  it  was  shown  that  there  is  always  at  least  one  solution  to  (1.1),  (1.2).  If 
A  “  1,  B”  0,  the  number  of  solutions  tends  to  “  as  e  ♦  0:  in  addition  to  the 
solution  y  «  -t,  for  each  positive  integer  j  and  e  =  e(j)  sufficiently  small,  there 
are  at  least  two  distinct  solutions  which  cross  y  -  -t  exactly  j  times  in  (-1,0).  The 
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limiting  behaviour  of  these  solutions  is  described  by  step  functions.  For  A  *  1,  B  =  0, 
necessary  and  sufficient  conditions  for  existence  of  solutions  were  derived  in  terms  of  the 
solution  to  an  associated  set  of  algebraic  equations.  The  main  techniques  involved  a 
theorem  of  Kabinowitz  14]  and  some  asymptotic  estimates  on  the  solutions  of  (1.1).  The 
restriction  B  =  0  was  necessary  mainly  in  order  to  use  Rabinowitz's  result. 

In  this  paper,  we  extend  the  results  of  [3]  by  removing  the  restriction  B  =  0.  In 
order  to  do  this,  we  use  completely  different  methods.  [3]  was  based  on  an  abstract 
bifurcation  analysis  which,  in  turn,  was  based  on  degree  theory;  this  paper  gives  a 
complete  analysis  based  entirely  on  a-priori  estimates  and  the  "shooting"  method.  In 
Section  2,  the  asymptotic  behaviour  of  solutions  is  analyzed  for  the  regions  |y|  <  |t| 
and  |y|  >  |t|.  It  is  shown  that  the  asymptotic  limits  of  the  crossing  points  of  y(t,£) 
with  y  ”  ±  t  can  be  explicitly  computed.  Section  3  is  used  to  derive  an  associated 
algebraic  system,  satisfied  (in  the  limit)  by  the  intersection  of  solutions  with  y  »  ±t; 
it  also  proves  the  existence  of  a  unique  solution  to  this  algebraic  system.  In  Section  4, 
we  return  to  the  question:  for  given  A  and  B,  how  many  solutions  are  there  to  (1.1), 
(1.2)?  The  answer  is  given  implicitly  in  terms  of  the  solutions  of  the  algebraic  system. 
For  most  A,  B  and  c  sufficiently  small  the  number  of  solutions  to  (1.1),  (1.2)  is 
finite  and  can  (with  some  work)  be  explicitly  calculated.  This  number  tends  to  “  as 
A  *  1,  B  *  0,  and  t  *  0.  If  A,  B,  the  number  of  "crossings"  j,  and  the  behaviour 

near  t  ■  -1  and  t  «  0  is  specified,  the  analysis  of  Section  4  yields  existence  and 

■uniqueness"  theorems  -  provided  e  is  sufficiently  small.  The  methods  of  analysis  enables 
one  to  get  a  complete  picture  of  the  variation  in  the  solutions  as  A  and  B  are  changed. 

Since  -y(t,e)  is  a  solution  of  (1.1)  whenever  y(t,e)  is,  we  may  restrict  ourselves 

to  the  case  B  <  A.  (For  a  ■  B,  y  =  A  Is  the  only  solution.)  The  case  B  <  A  <  0  is 

relatively  easy  to  handle,  and  will  be  discussed  in  an  appendix.  Thus,  our  main  effort  is 
concerned  with  the  case  B  <  A,  A  >  0.  Our  major  tools  are  those  developed  in  Section  2, 
using  extensions  of  the  asymptotic  estimates  of  [3)  and  backwards  shooting  methods  starting 
at  t  ■  0.  (The  backward  integration  is  used  for  a  technical  reason:  in  this  direction, 
it  can  be  proved  that  solutions  continue  to  exist  for  all  t.) 
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2.  Local  Results 


The  first  proposition  is  an  algebraic  result  which  will  be  useful  in  establishing 
facts  about  (1.1).  The  constant  t  satisfies  -1  <  t  <  0. 

Proposition  2.1:  Let  z(w)  =  |  w^  ■  t2w.  Then 

i)  z(w)  is  an  increasing  function  of  w  for  w  >  -t. 

ii)  If  y  >  -t  and  c  >  0,  there  is  a  unique  solution  w  >  -t  to 
z(w)  -  z(y)  +  c. 

Proof ;  i)  z(w)  is  a  cubic  with  critical  points  at  w  =  ±t.  The  point  w  *  -t  >  0  is 
local  minimum,  so  z(w)  is  increasing  for  w  >  -t. 

ii)  By  i),  y  >  -t  =*>  z(y)  >  z(-t).  Hence  for  c  >  0,  the  horizontal  line 
z  =  z(y)  +  c  is  aoove  the  local  minimum,  and  so  it  intersects  the  graph  of  z(w)  at  a 
unique  value  w  >  -t.  (See  Figure  2.1.)  I 


The  cubic  j  ■  j»  -  t2w,  t  <  0.  The  horizontal  line  is  z  -  z(y)  -  K,  y 
and  K  given.  The  desired  solution  to  z(w)  «  z(y)  -  K  is  the 
intersection  of  the  horizontal  line  with  the  darker  portion  of  the  cubic, 
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Let  y(t,e)  be  a  solution  to  (1.1). 

The  next  two  lemmas  describe  the  limiting  forms  of  the  trajectories  when  the  initial 
conditions  are  in  regions  I  and  11  respectively.  (See  Figure  2.2.) 


The  regions  (I):  y  >  |t|,  t  <  0  and  (II):  |y|  <  |t|,  t  <  0.  In 

region  I,  solutions  (integrated  backward)  have  derivatives  whose  absolute 
value  decreases  exponentially;  in  region  II,  the  derivatives  increase 
at  an  exponential  rate. 

-  -  -  _  i 

Lemma  2.1.  Suppose  that  lim  y(t,e)  -  y,  with  y  >  -t.  Suppose  also  that 
_  C*0 

lim  -  ey'(t,e)  =  X  >  0.  Then 
£♦0  _ 

i)  lim  y(t,E)  “  a  for  -a  <  t  <  t,  where  a  is  the  solution  >  -t  to 
e*0 

1—33—2— 

J  (y  -  a  )  -  t  (y  -  a)  »  -K  ,  (2.1) 

and 

-  r3  ,3 

ii)  lim  e  fn|y*(-a,E)|  -  -at  +  ^ - —  .  (2.2) 

6*0  J  J 

Proof:  we  first  show  that  y(t,E)  tends  to  a  limit  for  y  <  t  <  t,  some  yet.  (See 
Figure  2.3.)  Thi3  follows  from  the  direct  integration  of  (1.1):  We  get: 

.  t  _ 

y'(t,E)  »  [exp(-  jf  ( y  -  s  ) ds) ]y’ ( t,E )  .  (2.3) 

We  claim  that,  for  any  6  sufficiently  small,  lim  y’ ( t  -  5,e)  »  0.  For  y  >  -t,  let 

£♦0 
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0<i<y+t.  Then 


/  (y2  -  s2)ds  >  /  (y2  -  (t  -  5)2)ds  •=  o[y2  -  (t  -  6 ) 2 J  =  6  . 

t-5  t-6 


The  limiting  form  of  a  trajectory  with  initial  (t,y)  in  region  I,  and 

lim  -  ey' ( t,e )  =  K  >  0. 

E+0 


Thus  ly'(t  -  5,e)|  <  e  ly*  (t,e)  | ,  so  y’(t  -  i,E)  is  exponentially  small.  Now  suppose 
y  ■  t.  We  shall  show  that  we  can  replace  the  initial  conditions  t,  y  by  tE,  y(tE) 

(i.e.  a  point  on  the  same  solution)  with  some  tE  ♦  t  as  e  ♦  0  and  lim  y(tE)  >  y;  the 


previous  argument  will  then  hold.  First,  we  may  assume  that 


_  e*°  - 

t  <t,t  ♦  t  as  e  ♦  0 


such  that  y'(tE,e)  «  -K/2e  and  |y‘ |  >  K/2e  for  t  (tE,t] .  (If  not,  there  is  an 
interval  [t  -  a,t]  for  which  |y'(t,£)|  >  K/2 E;  it  is  then  clear  that  t,y(t)  may  be 

chosen  along  the  trajectory  such  that  lim  t  =  t,  lim  y  >  y.)  Next,  we  see  from  (2.3)  that 

e+o  e»0 

r  *  y' (tE,e)/y' (t,e)  =*  exp{-  ~  /  (y2  -  s2)ds}  . 

tE 

2  2  C  “ 

Either  y  -  s  is  unbounded  for  t  (t,t)  as  E  0,  in  which  case  we  are  through  (as 

above),  or  else  y2  -  s2  <  0(1).  In  the  latter  case,  since 


/  (y  -  s  )ds  "  c  In2 
t 
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we  must  have 


t|  >  0(C).  Finally 


e 

E  t  K 

y(t  ,e)  *  y(t,e)  +  /  y*(t,c)dt  >  y(t,e)  +  —  |t  -  t|  . 
t 

Hence,  lim  y(te,£>  >  y. 
e+0 

2  2 

As  long  as  y  -  t  >  0,  |y'(t)|  decreases  (integrating  backwards),  so  y'(t,£) 
remains  small  and  y  approaches  a  constant.  Since  this  is  true  for  any  small  6,  y  must 
approach  a  constant  for  y  <  t  <  t,  some  y.  In  order  to  see  what  that  constant  is,  we 
integrate  (1.1)  from  t=t  to  t  *  t  -  4.  Then 


e[y'(t  -  4 ,e )  -  y'(t,£)] 


u  2  2 
/  y'(s)(y  -  s  )ds  . 

t-S 


The  R.H.S.  of  (2.4)  may  be  rewritten: 


(2.4) 


t  ,  t  t , 

_  +  s  y(s,e)|_  -2  /  sy(s,e)ds  .  (2.5) 

t-4  Jt-4  -_6 


By  hypothesis,  and  by  the  previous  calculation,  the  limit  of  the  l.H.S.  is  K.  The  R.H.S. 
approaches 

3  -3  _ 

- — !j-*-  +  t  (y  -  a]  +  0(4)  .  (2.6) 

Since  (2.5)  is  valid  for  all  4  sufficiently  small,  we  conclude  that  a  must  satisfy 
(2.1).  (Note  that  a  is  the  solution  to  z(w)  =  z(y)  +  K,  which  exists  by  Proposition 
2*1.)  Furthermore,  it  is  now  clear  that  Y  may  be  taken  to  be  -a. 

ii)  To  establish  (2.2),  we  again  integrate  (1.1)  from  t  «  t  to  t  «  -a;  this  time, 
we  first  divide  by  y' .  We  get 


e  lnfy'(-a,e)|  -  e  £n|y'(t,e)| 


2  2 

"  /  (y  -  sz)ds  . 


(2.7) 


Since  y(t,e)  ♦  a  and  e  lnly'(t,t)!  *  0,  the  limiting  form  of  (2.7)  is  (2.2).  I 
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Ler.-ia  2.2.  Suppose  that  lim  y(t,E)  =  y,  with  |v|  <  |t|.  Suppose  also  that 

_  e+o 

lira  -  E  in|y'(t,E)|  =  L  >  0.  Then 

£♦0  _  _ 

II  lim  y(t,c]  =  y  for  -b  <  t  <  t,  where  b  is  the  unique  solution  >  lyl  of 

C+0 


-2-  t  -2,  b 

yt-  —  +  y  b  -  —  =-L  .  (2.8) 

ii)  3t£  *  -b  as  e  ♦  0  such  that  y  ( t^_ ,  £ )  =  -t£. 


iii) 


lim  e y '  ( t  , e )  =  — — — - b2(b  -  y) 

£♦0 


(2.9) 


Proof;  i)  Once  again,  we  integrate  (1.1)  to  get  (2.3).  Since,  by  hypothesis,  y’ (t,E)  is 
exponentially  small,  y* (t,e)  is  exponentially  small  for  all  t  >  -b,  where  b  satisfies 

t 

/  <y  -  s  )ds  *  -L  .  (2.10) 

-b 


(See  Figure  2.4.)  We  can  use  (2.10)  to  solve  for  b:  Equation  (2.10),  when  integrated, 
yields  (2.8).  Note  that  (2.8)  may  be  written  as 


J - y2b  *=  - y2(-t)  +  L  .  (2.11) 

By  the  same  arguments  as  in  Proposition  2.1,  with  y  and  -t  interchanged,  there  is  a 
unique  solution  b  to  (2,11)  for  b  >  y. 

ii)  For  any  S  >  0  sufficiently  small,  (2.3)  shows  that  y' (-b  -  6 , c )  either  grows 
exponentially  large  in  e,  or  else  3t£,  with  -b  -  5  <  t£  <  -b,  such  that 

y(t£,e)  “  -t£.  (See  Figure  2.3.)  Since  the  exponential  growth  also  implies  the  existence 
of  such  a  t£,  and  since  this  is  true  for  all  6  >  0,  assertion  ii)  is  proved. 

iii)  The  growth  of  y  implies  that  y' (t,e)  must  stop  being  exponentially  small  in 


c,  i.e.  3  t£  >  t£,  with  t£  *  b  as  e  +  0,  such  that  y' (t£,e)  ■  -e.  Note  that, 

since  y'(t,E)  ♦  0  for  t  <  t  <  t,  lim  y(t,c)  ■  y  for  t  <  t  <  t.  To  establish  (2.9), 

E  £>0  C 

we  now  use  equations  (2.4)  and  (2.S),  with  t  -  5  and  t  replaced  by  t  and  t£ 
respectively.  ) 


-7- 


y 


Figures  2.4a,  2.4b 

The  limiting  form  of  trajectories  with  initial  t,  y  in  region  IX,  and 
lim  -  e  Zn(y*(t,e)|  =  L  >  0.  Figure  2.4a  is  for  y  >  0;  Figure  2.4b 
il  ?or  y  <  0. 


The  last  lemma  of  this  section  describes  the  initial  segment  of  the  trajectory  of 
(1.1)  if  y(0,E)  -  b  <  a. 


Lemma  2.3.  Suppose  lim  y(0,e)  =  B  <  0,  and  lim  -  ey*(0,e)  =»  K  >  0.  Then 
e-*0  e-*0 

i)  lim  y(t,e)  =  a  for  -|a|  <  t  <  0  ,  where  a  is  the  real  solution  to 
e-»0 

i  a  i 

K  , 


j  (a3  -  B3) 


(2.12) 


and 


ii) 


lim  e  ln|y* { —  I  a | ,£) !  =  -  \  |a|3  . 
e+0  J 


(2.13) 


Proof ;  This  lemma  is  essentially  the  same  as  Lemma  2.1,  with  t  =  0i  there  is  an  added 
difficulty  that  the  trajectory  may  cross  through  the  region  |y|  <  |t|,  so  y*  may  not 
change  in  a  monotone  fashion.  (See  Figure  2.5.)  Hence  we  cannot  immediately  use  the 
previous  argument  that  y(t,e)  approaches  a  limit  for  t  near  t  «  0. 


-8- 


I 


The  limiting  form  of  trajectories  for  initial  condition  y(0,s)  =  B  <  0, 
and  lim  -  ey'(0,e)  =  K  >  0.  In  Figure  2.5a,  K  is  small  enough  that 
lim  y(?,e)  <  0  for  t  near  0.  If  the  limiting  form  is  as  in  Figure  2.5b, 
for  e  >  0  the  trajectories  pass  through  |y|  <  ]t). 

Let  a  be  defined  as  in  (2.12),  and  let  5  >  0  be  arbitrarily  small,  we  first  claim 
that  y(t,e)  cannot  remain  bounded  away  from  v  «  a  for  -5  <•  t  <  0  and  all  e,  6 

sufficiently  small.  For,  integrating  (1.1)  from  t  =  0  to  t  »  t»  (-5  <  t#  <  0 )  ,  we  get 

0  2  2 

C ly * (t4,e)  -  y'(0t£)]  =  -  /  y’(s)(y  -  s  Ids  .  (2.14) 

By  hypothesis,  -lim  ey'(0,E)  =  K  >  0.  If  y  stays  bounded  away  from  a,  the  R.H.S.  of 
e->0  _  _ 

(2.14)  will  be  strictly  less  than  K  +  0(5).  Thus  !y'(t#,e)|  >  K/e  for  some  K  >  0  and 

6  sufficiently  small.  Since  this  is  true  of  all  t»  such  that  -5  <  t#  <  0,  y  changes 

*>y  0(1/e),  and  cannot  stay  bounded  away  from  a  for  any  finite  value  of  t;  this 

contradicts  the  assumption  that  a  is  an  upper  bound  for  y. 

Also,  y  cannot  become  strictly  larger  than  a,  i.e.  bounded  away  from  a  from 
above  for  any  t  <  0  with  t  sufficiently  small.  For,  if  y  hits  y  =  a,  let  t^(e) 
be  the  value  of  t  when  y  «  a.  Then  from  (2.12)  and  (2.14)  we  have  that 
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lim  ty'(t  ,£)  =  0.  For  all  t  such  that  -|a|  <  t  <  0,  |y‘ |  decreases  (as  (1.1)  is 
integrated  backwards),  and,  as  in  Lemma  2.1,  y(t,£>  approaches  a  as  £  ♦  0. 

ii)  Same  as  the  proof  in  Lemma  2.1.  I 

Corollary  2.1.  Suppose  lim  -  £y* (0,£)  =  K  >  0,  and  lim  j  lim  y(S,£)  =0  or  “.  Then 

£♦0  &+0  t* 0 

y(t,£)  exists  for  all  t  <  0.  Furthermore,  if  the  solution  hits  y  *  -t  at  t  <  0,  and 

lim  £  in  y'(t,e)  *  0  at  the  crossing  (if  that  crossing  is  horizontal)  or 
e+0 

lim  ey' (t,e)  *  0  (if  the  crossing  is  vertical),  then  y(t,£)  crosses  y  »  -t  an 
£■*•0 

infinite  number  of  times. 

Proof :  If  B  >  0,  we  apply  Lemmas  2.1  and  2.2  alternately,  starting  with  Lemma  2.1  (since 

y  >  t  =  0).  These  lemmas  explicitly  calculate  the  limiting  values  of  the  successive 

crossing  points.  (If  K  =  0,  there  may  be  only  one  crossing  point,  e.g.  if 

y‘(0,£)  =  0  Ve . )  If  B  =  0,  the  hypothesis  on  lim  4-  lim  y(5,e)  forces  .he  behaviour 

£•*•0  e>0 

near  t  *=  0  to  be  either  horizontal  or  vertical;  we  can  then  apply  Lemmas  2.1  and  2.2. 

If  B  <  0,  the.  above  hypothesis  again  rules  out  solutions  y(t,e)  for  which  y  *  -t  is 
the  limiting  solution  for  -1  <  t  <  0.  We  can  start  with  Lemma  2.3  and  then  apply  Lemmas 
2.2  and  2.1  alternately.  The  uniqueness  of  the  point  B(e)  is  clear  from  thse  lemmas.  | 
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3.  An  Aluehraic  Svrtem 


We  have  seen  that,  in  the  Unit  as  c  ♦  0,  almost  all  of  the  solutions  to  (1.1)  are  a 

succession  of  horizontal  and  vertical  segments;  the  length  of  each  segment  can  be 

explicitly  computed  given  lim  y(0,£)  and  lim  cy'(0,c). 

e*0  £+0 

We  shall  consider  a  system  of  algebraic  equations  which  are  satisfied  in  the  limit  by 
the  successive  points  t^  at  which  y(t,e)  crosses  the  line  y  =  -t.  Much  of  what  is  in 
this  section  was  already  done  in  [3)  under  the  condition  that  B  =  0.  Here  we  show  that 
those  results  still  hold  for  3*0.  The  first  lemma  and  theorem  apply  to  a  solution 


y(t,£)  if  lim  y(t,£)  >  0  for  t.  <  0. 

£•*•0 

Lemma  3.1.  Let  t.  .(£)  >  t.(£)  >  t.  .(£)  be  three  successive  points  at  which  a  solution 

-  1-1  i  1+1 

y(t,£)  crosses  y  =  -t.  If  t^._i ,  t^,  t^+^,  denote  the  limits  of  these  points  as 
£  +  0,  then 


.2 

'I 


J  (ti-i  +  ViVi  +  Vi’ 


(3.1) 


Proof;  There  are  two  possible  cases,  as  shown  in  Figures  3.1a,  3.1b.  Each  case  can  be 
done  by  the  methods  of  Section  2.  For  the  case  of  3.1a,  let  6  >  0  be  any  sufficiently 
small  number,  and  use  (2.4)  with  t  -  5  and  t  replaced  by  t^  -  6  and  t^  +  5.  Since 
y(t^  +  6 ,  £ )  +  ”ti_1  and  yt^  -  i,E)  *  _ti+ir  using  (2.5)  and  letting  6  *  0  we  get 
(3.1). 

For  the  case  of  3.1b,  we  use  (2.7),  with  -a  and  t  replaced  by  ti+1  and 
By  previous  computations,  the  L.H.S.  tends  to  zero  as  E  +  0.  Also,  for  ti+1  <  t  < 
y(t,£)  +  t^.  Evaluating  the  R.H.S.  and  setting  it  equal  to  zero,  we  again  get  (3.1).  I 
We  shall  now  prove  an  existence  and  uniqueness  theorem  for  (3.1),  regarded  as  a  system 
of  equations  for  a  set  of  points  {t^} .  We  shall  denote  by  tg  and  t^  the  right  and 
left  endpoints  of  the  set. 

Theorem  3.1.  There  is  a  unique  solution  0  >  tQ  >  t^  >  •••  >  t^  to  (3.1)  with 


given 


(3.2) 


Figure  3.1 

Successive  crossings  of  y  =  -t  by  limits  of  solutions  to  (1.1)  as 
C  ♦  0. 


Proof;  Equation  (3.1)  gives  rise  to  a  mapping  of  finite  sequences  of  numbers  between  t- 


and  0,  namely;  For  any  sequence  »  V  <  <  •••  <  xQ  »  tQ,  let 

VT>  --{T  tTli  +  V1V1  +  Ti+i^}  1/2  (3*3) 

for  i  *»  1,2, ...,j,  and  FQ  =  tQ,  F^  =  t  j .  It  is  clear  that  any  fixed  point  of  this  map 
is  a  solution  to  (3.1)  and  (3.2).  We  shall  show  that  {f^}  has  a  unique  fixed  point. 
Consider  two  such  sequences  T  =  {tJ  and  0  “  {o^}.  Then 

Fj(a,  -  fJ(T,  -  1  ((02i+1  -  T2i+1)  +  (°i+iai_l  -  'it1Vl'  +  <°i-1  -  *!-,»  (3-4) 


for  i 


1.  Let 


=  o .  -  T,  W .  =  0,  +■  T..  Also,  write  W.  ,  for 

i  i  i  r  i  i  r  i 


(3.5) 


F^to)  +  TA"Z)  and  F^ta)  -  F^(  T ) .  Equation  (3.4)  may  then  be  written  as 

Vi  *  J  (wi+i  +  1  wi-i^i+i  +  7  (wi-i  +  7  Vi’Vi  ' 

i  •  l,...,j  -  1.  AlSO,  ^  "  0,  •  0. 

The  coefficients  of  in  (3.5)  are  not  constant,  but  we  will  get  bounds  on  them 

which  insure  that  {f^}  maps  any  two  sequences  closer  to  one  another.  We  divide  (3.5)  by 


W^,  and  write  the  resulting  system  as 


14 


(3.6) 


Notice  that  the  matrix  L  has  non-zero  entries  only  just  above  and  just  below  the 


diagonal.  We  first  obtain  a  bound  on  those  coefficients,  and  then  show  this  implies  the 


desired  result.  Let  a^  denote  the  entry  in  the  i*"*1  row,  (i  +  1)t^'  column;  is  in 


the  it^'  row,  (i 


1  ) th  column. 


Lenma  3.2. 


Uil  +  1  bi  |  <  1. 


Proof ;  a .  =  — —  [W .  .  +  —  W .  .  1  :  b. 

-  l  1-1  2  1+1  i 

3W. 

An  elementary  calculation  shows  that 


3W, 


Vi  +  2  Vi1 


a.  +  b. 
i  i 


Now 


Hence 


2W. 

1 


Vi  +  wi-i] 


Both  a^  and  b^  are  >  0. 


(3.7) 


0  <  h  Vi  -  Ti-i>2 


,1  In  2  ,  1  2n  , 1  In  2 

(j  -  4  Vi  +  <3  -  4  VVl  +  ^3  "  4  Vi 


j  ‘v  *  Wi ♦  t2-i}  >  i  {Vi +  Wi +  Tli} 


This  mav  be  rewritten  as 


lFi{T)l  >  2  Vi  +  Vi1 


(3.8) 


Inequality  (3.8)  also  applies  to  o,  and  hence 


'V  >  2  'Vi +  Vi1  * 

It  follows  that  the  R.I1.S.  of  (3.7)  is  <  1 .  I 

We  return  to  Theorem  3.1.  L  is  a  tri-diagonal  matrix  whose  entries  satisfy 
|a^|  +  IbjJ  <  1.  It  follows  (6]  that  if  |X|  >  1,  then  L  -  XI  is  nonsingular,  so  all 
the  eigenvalues  of  L  are  <1  in  absolute  value.  From  this  it  can  be  shown  that  {f^} 
has  a  unique  fixed  point.  In  fact,  a  direct  iterative  scheme  is  convergent  to  this  unique 
fixed  point,  | 

We  now  deal  with  solutions  which  cross  y  m  +t  as  well,  i.e.  when  B  <  0.  We  have 

seen  that  there  is  a  unique  turning  point  0(E),  and  we  shall  be  concerned  now  with 

solutions  for  which  8  =  lim  8(e)  <  0.  Lemma  2.3  implies  that  for  t  <  0  sufficiently 

C*0 

near  t  »  0,  lim  y(t,e)  ■  8. 

£♦0 

Lemma  3.3.  Suppose  lim  y(t,e)  -  6  <  0  for  8  <  t  <  0.  Let  t1 (e ),..., t. (e) ,.. .  denote 
€♦0  3 
the  crossings  of  y(t,£)  by  y  -  -t,  and  t, ,t2, • • • , t^, . . .  their  limiting  values.  Then 


i) 


(3.9) 


3.2a  3.2b 

Figures  3.2a,  3.2b 

Crossings  of  y  =  +t  and  y  -  -t  by  a  limit  of  solutions  if  B  <  0  and 

2  2 

8  <  0.  For  Figure  3.2a,  8  and  t^  are  related  by  38  =  t^.  For 

Figure  3.2b,  t^  >  3&2 • 
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I 


Theorem  3.2.  1)  Let  t^  be  given.  Then  there  is  a  unique  solution 
0  >  8  >  t1  >  t2  >  •••  >  t  to  (3.1),  i  =  =  1,  and  (3.9),  (3.10). 

ii)  Let  tj  and  8  be  given.  Then  there  is  a  unique  solution  0>8>t,  >»**t 
to  (3.1),  i  =  2,...,j  -  1  and  (3.9). 

Proof :  The  proof  is  almost  identical  to  that  of  Theorem  3.1,  so  we  shall  just  indicate  the 
differences. 

i)  Let  T  denote  a  sequence  8>T,  >t2>*»»>t^,  0  a  sequence 
T  >  o,  >  >  By  F^d),  i  >  2  is  as  before,  but 

VT)  =  ^3  ‘  BT2  +  T2^ 


Fo(T)  VTi 


W.,  $  ,  W.  are  as  before.  Then  (3.5)  still  holds  for  i  >  2.  For  i  «  0,1,  we  have 


Vo  "  3  Vi  ' 


Vi  -?  (w2  -i  V*2  +  i  tw« 


(3.17) 

"0  "  2  V90  *  t3,12) 

To  show  that  la^J  +  ib^l  <  1  for  all  i,  it  suffices  to  show  this  for  i  *  1  and 
i  *■  0.  For  i  =  0,  bj^  »  0  and  la^l  ■  w^/3wg  ”  1/^3  <  1.  For  i  =  1, 


la, |  +  |b, | 


f  (|t2i  +  io2i  -  : e i  -  iy i ) 

|f,(t) |  +  |F,(o)i 


As  in  the  previous  theorem,  it  can  be  shown  that  |F,(t)|  >  |t  -  61  and  similarly  for 
F, (o) ,  which  implies  that  1  |  +  |b^|  <  1. 

ii)  Let 

F,d,  --{I[62  -  6t  +  x2]}  . 

Fi<T) ,  i  >  1  as  before  . 

The  relevant  equations  are  (3.5),  i  A  2,  (3.12).  The  proof  that  |a,|  <1  is  done  by 

the  usual  methods.  I 
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4  •  Existence  an»:  Vnicneni";::  of  Solutions 

The  solution  to  the  boundary  value  problem  (1.1),  (1.2)  is  not  unique,  as  shown  by 

examples  in  [1],  [3].  However,  as  we  shall  see,  any  two  solutions  are  qualitatively 

different:  they  can  be  distinguished  by  the  number  of  points  at  which  the  solutions  cross 

y  “  It,  and  the  behavior  of  the  solutions  at  the  boundary. 

From  the  estimates  of  Section  2,  it  follows  that  the  behavior  near  each  boundary  is 

characterized  by  lim  ey‘(t,c),  where  t  =  -1  or  0;  if  this  limit  is  non-zero,  the  initial 
e*0 

segment  of  the  solution  at  that  boundary  approaches  a  vertical  line  whose  length  goes  to 

zero  with  this  limit.  If  lim  y'{t,e)  =  0,  the  initial  segment  is  (in  the  limit)  horizontal. 

e*0 

We  shall  divide  the  analysis  into  four  cases,  depending  on  whether  the  behavior  at 
t  «  -1  and  t  =  0  is  vertical  or  horizontal.  In  each  case  we  aslc  the  question:  for 
given  A  and  B,  how  many  solutions  are  there  to  (1.1),  (1.2)  in  the  limit  as  e  *  0? 

The  answer  is  given  implicitly  in  terms  of  the  solutions  to  (3.1),  but  in  a  form  such  that 
it  could  be  explicitly  computed.  The  analysis  provides  a  complete  picture  of  how  the 
solutions  change  as  A  and  B  are  varied,  and  also  how  the  boundary  value  A  changes 
as  B  and  y' (0,e)  are  varied. 

We  first  discuss  solutions  which,  in  the  limit  as  e  *  0,  have  a  vertical  segment 

(shock  layer)  at  both  t  »  -1  and  t  *  0.  (See  Figure  4.1.)  We  also  assume  that,  if 

B  <  0,  then  8  «  0.  For  each  odd  j  we  consider  the  solution  to  (3.1)  with  tQ  -  0 

and  tj+i  ”  “I*  (For  this  bundary  behaviour,  the  number  of  strictly  interior  crossing 

points  must  be  odd.)  Let  A j,  j  >  1  be  the  open  interval  (-t^,  where 

denotes  the  itf'  point  of  the  solution  to  (3.1)  having  j  +  2  points  t^  <  •••  <  *  1  . 

(If  the  j  is  understood,  we  may  sometimes  suppress  the  superscript.)  Similarly,  let  3 j 

be  the  open  interval  (-”,-tj).  By  "strictly  interior"  points  of  intersection  with 

y  »  -t,  ve  mean  those  t^te)  for  which  -1  <  lim  t?(e)  <  0. 

1  e+0 

Remark :  For  fixed  c  small,  the  actual  number  of  interior  points  of  intersection  may  be 
greater  than  the  limiting  number  of  "strictly  interior"  such  points.  That  is,  we  Blight 


have  tQ(e)  <  0,  ( c )  ♦  0  as  c  *  0  or  t^+1(e)  >  -1,  lim  t^+1(e) 


The  above 


method  of  counting  was  adopted  because  it  eliminates  the  division  into  further  special 


-16- 


cases  that  would  be  needed  if  those  other  points  of  crossing  were  counted  as  well,  end 
simplifies  the  description  of  the  collection  of  all  solutions. 

Theorem  4.1.  i)  There  is  a  solution  to  (1.1),  (1.2)  having  a  boundary  layer  at  each  end, 

j  strictly  interior  points  of  intersection  with  y  =  -t  and  no  strictly  interior  points 

of  intersection  with  y  =  t  if  A  e  A^,  B  e  Bj  and  c  is  sufficiently  small.  Such  a 

solution  is  unique.  Its  limiting  behavior  as  E  ♦  0  is  a  step  function:  a  series  of 

horizontal  and  vertical  lines  crossing  y  =  -t  at  the  points  t?.  A  necessary  condition 

for  such  a  solution  is  A  £  c£  .4  ,  B  £  cl  3  ,  where  "cl"  denotes  "closure  of". 

1  3 


A  solution  to  (1.1)  with  boundary  layers  at  both  ends  and  three  interior 
turning  points. 

11)  The  {/O  and  £3^}  are  nested  intervals,  i.e.  ^  O  D  D  •  •  • ; 

fl  D  3  3  3  3  •••. 

1  3  5 

Hi)  For  each  A  *  1,  (resp.  B  >  0),  there  are  finitely  many  j  for  which  A  e  A ^ 
(resp.  Be  S  ^ ) . 

Proof;  i)  The  necessity  of  the  condition  U  cl  A  Be  cl  B  ^  follows  immediately  from 
the  results  of  Sections  2  and  3.  For,  by  Section  2,  the  limiting  form  of  a  solution  is  a 
sequence  of  horizontal  and  vertical  lines,  and,  by  Section  3,  the  intersections  of  these 
lines  with  y  «  -t  are  the  points  t^.  The  hypothesis  A  e  cl  4.,  s  «  cl  3^  requires  that  A 
and  B  lie  on  the  closure  of  the  vertical  segments  at  t  •  -1  and  t  *  0  respectively,  for  the 


limiting  solution  to  (1.1)  corresponding  to  the  solution  to  (3.1)  with  t  +j  =  -1#  tg  =  0; 
if  this  hypothesis  is  not  satisfied,  the  limiting  solution  cannot  be  the  correct  one. 

To  establish  the  existence  and  uniqueness  of  solutions,  consider  the  solutions 
y(t,£)  to  (1.1)  satisfying 


y(0 ,£)  --  b,  y'  (0,0  =  -k/e  . 


(4.1) 


By  the  estimates  of  Section  2,  there  is  a  unique  lO  *  0  such  that  the  solution  with 

K  *  Kj  has  the  desired  limiting  value  as  c  •*  0.  (K^  is  chosen  so  that 

lim  y(t,£)  =  -t3  for  t  near  0;  this  uses  B  f  2,.  The  rest  of  the  limiting  behavior 

€♦0 

follows  automatically  from  the  results  of  Section  3.)  We  shall  show  that  if  A  t  A ^  there 

is  a  unique  K^(£)  such  that  K^(c)  *  iO  as  e  +  0 ,  and  the  solution  through 

y(0,O  »  B  with  slope  y'(0,£)  *»  -K^(e)/£  satisfies  y(-1,E)  =  A.  (This  "local" 
uniqueness  suffices  to  insure  uniqueness  of  solutions  with  boundary  layers  at  both  ends 
and  j  -  1  intersections  with  y  a  -t,  since,  by  Section  3,  any  such  solution  must  have 
the  limiting  behavior  described  in  the  theorem.) 

Consider  the  solutions  to  (1.1),  (4.1),  with  |K  -  K^|  <  5,  5  arbitrarily  small.  The 
next  two  lemmas  are  monotonicity  results.  The  first  one  concerns  the  limiting  values  of 
the  solution,  and  could  be  phrased  as  a  statement  about  equations  (3.1).  (See  Lemma  4.2.) 


Lemma  4.1.  Let  t^(£,K)  be  the  value  of  t  where  the  solution  y(t,£,K)  to  (1.1),  (4.1) 

hits  y  »  -t  for  the  3th  time.  (By  definition,  lim  t  (£,K)  «  -1.)  Then  for  5 

£♦0  3 

sufficiently  small. 


dt 


lim  — ' i  (E.K)  <  0  , 

e+o  dK 


(4.2) 


so  lim  t.(e,K)  is  a  monotone  decreasing  function  of  K. 
e*0  3 

Proof:  The  monotouicity  follows  from  the  estimates  in  Section  2.  For,  using  the  cubic  of 

Proposition  2.1,  it  can  be  seen  from  Lemma  2.1  that  increasing  K  or  increasing  B  has 

the  effect  of  increasing  -lira  t  (e,K)  and  -lim  £  ln|y' (t  ,£ ,K) | , 

£♦0  c+o 


in  a  transversal  way: 


Similarly,  from  Lemma  2.2,  increasing  lira  |t  (£,K)|,  lin  |c  in | y 1 ( t  , £ , K) | |  increases 

£-«•  o  1  e-*0  1 

lim  |t  (£,X)|  and  lim  £  I y ' ( t  ,£ ,K) | .  Continuing,  alternating  Lemmas  2.1  and  2.2,  we 
£■*0  2  £-*0  ~ 
find  that  (4.2)  holds.  I 

We  now  return  to  the  proof  of  Theorem  4.1.  Lemma  4.1  implies  that  for  K  <  (resp. 

K>  Kj),  t^+1(£,K)  >  t^+1(c,K^)  (resp.  t^+1(£,K)  <  t^te.K^))  for  all  e  sufficiently 
small.  Thus,  for  £,  6  small,  y(-1,E,K)  is  arbitrarily  close  to  the  upper  limit  of  A ^ 
for  K  <  Kj  and  the  lower  limit  of  A ^  for  K  >  (see  Figure  4.2).  As  we  will  see, 

y(-1,e,K)  is  not  monotone  for  |K  -  K^l  <  6,  any  fixed  6,  for  all  e  sufficiently 
small.  However,  we  shall  show  that  there  is  a  smaller,  e-dependent  interval  in  K  such 
that  y(-1,e,K)  is  a  monotone  function  of  K  on  this  interval,  and  the  image  of  this 
interval  tends  to  as  e  ♦  0.  Furthermore,  as  K  passes  the  boundaries  of  this 

interval,  the  size  of  the  "vertical"  segment  at  t  =  -1  tends  to  zero,  so  the  limiting 
behavior  changes  from  vertical  to  horizontal.  This  suffices  to  prove  the  uniqueness  of  the 
solution  with  boundary  layers  at  both  ends,  and  j  (strictly)  interior  turning  points. 


The  limiting  form  of  nearby  solutions  to  (1.1)  with  the  same  value  of 
B  and  different  values  of  y' (0,e)  «  -K/e,  (j  “  1).  For  K  »  K^, 
the  limiting  solution  has  a  vertical  segment  at  t  -  -1 .  For  K  <  K ^ 
(resp.  K  >  K j )  the  limiting  solution  has  a  horizontal  segment  near 
the  upper  (resp.  lower)  endpoint  of  ,4^. 
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t 


By  Lemma  4.1,  for  each  e  sufficiently  snail,  t_.  +  i)(e,K) 

K  near  K^,  so  31  K^(£)  such  that  t  ( c  ,K^  ( £ ) )  =  -1.  By 

y(tj+1<e,K),e,K)  =  -t..+  1(e,K)  . 

Differentiating  (4.3)  with  respect  to  Kj  and  evaluating  at  K  ° 


is  a  monotone  function  of 
definition  of  t^+^(e,K), 
(4.3) 

“  K  (e) ,  we  get 


dt . 


4*  <-1,e,k  <£))  =  -[1  +  y‘  (-1  ,  £  ,K  (£ ) )  ]  -  (£,K^J 


+  1 


3k 


dK 


3 1 


(4.4) 


Now  y'  (t  ^ (£,K)  ,£,K)  is  <0  and  0(1/e).  Furthermore,  by  Lemma  4.1, 


dt^+1 (e,K)/dK  <  0.  Hence  (3/3K)y(-1  ,£  ,K^  ( £  ) )  <  0  and  0(1/£),  and  there  is  an 


interval  K_._(£)  <  K  <  K_.^ { e  }  such  that  ( 3/3K)y(-1 ,e ,K)  <  0. 

It  remains  to  be  shown  that  the  above  interval  may  be  chosen  so  that  the  image  under 


the  map  K  y(-1,£,K)  tends  to  A ^  as  e  ♦  0,  and  that  outside  this  interval,  the  limit 


behavior  near  t  =  -1  is  horizontal. 

To  see  this,  we  note  that  the  monotonicity  persists  as  long  as  ( o/3K)y( -1 ,£  ,K)  <  0. 
Also,  the  conditions  y  *  y(t,e,K),  (3/3K)y( t,e ,K)  =  0  define  the  envelope  of  the 
trajectories  y(t,£,K)  parameterized  by  K.  Thus,  monotonicity  persists  until  a  curve 
y(t,£,K)  touches  the  envelope  at  t  «  -1.  From  the  variational  equation  of  (1.1)  around 


y(t,£,K  ),  it  can  be  seen  that  ( 3/3K)y( t,£ ,K) I  <  0  for  1 1  —  1 1  <  0(/F).  Since 

jl  K-K^ 

the  width  of  the  boundary  layer  is  0(e),  it  follows  that  the  curves  y(t,e,K)  for 


nearby  K  intersect  near  the  top  and  bottom  of  A ^  for  e  small.  (Geometrically,  the 


envelope  of  the  limiting  functions  (as  e  ♦  0 )  is  a  pair  of  curves  through  the  "corners" 
separating  the  horizontal  and  vertical  segments.  For  £  >  0,  the  envelope  is  a  pair  of 


curves  converging  to  the  limiting  pair.)  Outside  the  interval  K^_(£)  <  K  <  K_.  +  (£),  and 


for  |K  -  K^ |  sufficiently  small,  the  size  of  the  "boundary  layer  jump" 


lim  lim  I y ( — 1 ,e,K(£) )  -  y(t,£,K(£})|  is  zero. 
t*-1  e*0 

Having  proven  i),  we  now  turn  to  the  proof  of  ii ) .  We  first  show  that 


3  it  suffices  to  prove  that  t^+2  >  t^. 

Lemma  4.2.  Suppose  t^  <  t^.  Let  {t^},  {t^}  be  the  solution  to  (3.1)  with  tg  “  -1, 


■  t^  or  t1,  respectively.  Then  for  each  i,  t^  <  t^. 
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Proof :  31  K  (resp.  K)  such  that  the  solution  to  (1.1)  satisfying  y(0,E)  -  0, 
y'(0,C)  ■  -c/K  (resp.  -c/K)  crosses  y  «=  -t,  in  the  limit  as  e  +  0 ,  at  t1  (resp. 
t1 ) ;  if  t  <  t  ,  then  K  >  K.  We  know  from  Lemma  4.1  that  the  successive  crossing 
points  are,  in  the  limit,  monotone  functions  of  K,  and  from  Section  3  that  these  crossing 
points  approach  solutions  of  (3.1).  I 

We  now  return  to  the  proof  of  ii).  It  follows  immediately  that  t^+^  >  t^  .  For  if 

not,  we  would  have  t^+]  >  t^+^,  for  some  j;  but,  by  hypothesis,  t^+]  =  -1  and  t?+^  >  -1. 

g+1  3+1  r  3+1  3+1 

A  similar  argument  shows  that  3  "j+2"  T^at  is,  the  point  -1  and  the  next 

nearest  points  t^+1  or  tj+2  determlrie  a  solution  to  (3.1).  As  in  Lemma  4.2, 

t"!+1  >  t^*i?  and  tj+’  <  t^+^,  so  4.  D  A.  .. 

3  3+2  3+2  3+4'  3  3+2 

In  order  to  prove  iii),  first  consider  B  >  0.  Since  t^  is  monotone  in  j,  it 

suffices  to  know  that  t^  +  0  as  j  ♦  since  B  will  then  eventually  be  outside  3j 

for  large  enough  j.  This,  in  turn  follows  by  considering  solutions  to  (1.1)  which,  in  the 

limit  as  e  +  0,  crosses  y  **  -t  for  the  first  time  at  t ^  •  It  is  then  clear  that  if 

t^  is  bounded  away  from  0,  we  could  not  have  the  required  “  -1.  (Indeed,  we 

j+  1 

would  have  t  .  +  -*•.) 

J+1 

A  similar  argument  holds  for  A  "  Is  i.e.  t^+1  +  ~1»  t^*^  *  -1»  so  A  .  1  must 


eventually  fall  outside  A ^  for  sufficiently  large  j. 


I 


Note  that  this  argument  does  not  rule  out  a  countable  number  of  "limiting  solutions" 
for  A  ”  1,  B  <  0.  However,  the  bound  on  e  is  necessary  to  insure  the  existence  of  a 
real  solution  with  j  strictly  interior  crossings  with  y  »  -t  may  go  to  zero  as  j  +  •». 

Next,  we  deal  with  solutions  (other  than  the  constant  solutions)  which,  in  the  limit 
as  c  ♦  0,  have  horizontal  segments  near  both  t  »  -1  and  t  ■  0. 

Theorem  4.2.  Suppose  B  >  0.  For  each  odd  j  >  3,  consider  the  solution  to  (3.1)  with 

A  (If  A  >  1),  and  t,  -  -B  (if  B  >  0)  or  tQ  =  -B 
(if  B  -  0).  (See  Figures  4.3a  and  4.3b.) 

i)  There  is  a  unique  solution  to  (1.1),  (1.2)  which  approaches  horizontal  segments 
at  t  «  0  and  t  «  -1,  and  which  has  3  strictly  interior  turning  points  if  the 
following  is  satisfied!  if  A  <  1,  then  t^+^  <  —1  *  if  A  >  1,  then  -1  <  tj  <  0.  Also 


tj  ■  -A  (if  A  <  1 )  or  t 
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(tj)^  >  3 A  necessary  condition 

2  2 
(tj)  >  3B  . 


for  such  a  solution  is  t .  .  <  -1  if  A  <  1  and 

3+1 


Figures  4.3a,  4.3b 

Solutions  to  (1.1)  with  j  =  3  and  "horizontal"  segments  near  t  *  -1 
and  t  *=  0.  Figure  4.3a  shows  A  <  1;  Figure  4.3b  has  A  >  1. 

.2  .2 

ii)  If  for  some  j,  (t^j  <  332,  then  (t^*2)  <  3B2.  If  A  <  1  and  tj+1  >  _1» 

then  t^*2  >  -1.  If  A  >  1  and  f  (-1,0),  then  t^2  t  (-1,0). 

iii)  If  A  *  1  or  B  >  0,  there  are  only  finitely  many  solutions  to  (1.1),  (1.2). 

Proof:  i)  We  first  establish  the  necessity  of  the  conditions  in  i) .  As  in  Theorem  4.1, 

the  conditions  come  from  the  requirement  that  the  limiting  solutions  intersect  y  =»  -t 
exactly  at  the  {t?},  where  {t:|}  is  the  solution  to  (3.1)  with  the  boundary  conditions 
given  in  i) . 

2  2 

The  condition  (t^)  >  3b  requires  that  there  be  no  negative  solution  to  the  i  =  1 
equation  of  (3.1)  (with  t1  »  -B,  t2  given  by  the  appropriate  solution  to  (3.1)). 
Equivalently,  there  is  no  "next  point"  to  the  right  of  t^  which  is  still  <  0.  (For  the 
limiting  solution  to  (1.1),  (1.2)  must  have  -B  as  the  first  point  of  intersection  with 
y  ■  -t  if  B  >  0 j  for  B  «  0  the  condition  is  vacuous.)  Similarly,  the  condition  on 
tj+j  (if  A  <  1)  demands  that  t  =  -A  be  the  last  (limiting)  point  of  intersection 
before  t  -  - 1 j  if  a?  1,  the  condition  on  t^  requires  that  t  »  -A  be  the  first 


point  of  intersection  for  t  <  -1 


Wc  now  90  to  the  existence  and  uniqueness  of  the  solutions.  Using  Lemmas  (2.1)  and 


(2.2)/  we  see  that  (t ^ )2  >  3B^  implies  3(  L.  *  0  such  that  the  solution  to  (1.1)  with 


y(0,e)  »  B,  y '  (0, e )  =  -e 


-L/e 


(4.4) 


and  L  =  Lj  has  a  limiting  behavior  whose  first  vertical  strip  is  through  t^.  (It 

follows  that  the  limiting  behavior  of  this  y(t,£,L_.)  is  a  series  of  horizontal  and 

vertical  strips  cutting  y  ■=  -t  through  all  the  {t^}.)  We  shall  show  that  for  each  e 

sufficiently  small,  2  J  o  (c)  such  that  L  (£)  L .  as  e  *►  0/  and  the  solution  through 

3-L  (£)/c  3  3 

y(0,e)  =  B,  y*  (0,0  =  -e  ^  satisfies  y(-1,e)  =  A.  As  in  Theorem  4.1,  we  need  a 

monotonicity  result. 

Lemma  4.3.  Let  t^(e,L)  be  the  value  of  t  where  the  solution  y(t,e,L)  satisfying 

*■  h 

(4.4)  hits  y  =  -t  for  the  i'~1  time.  (By  definition,  lim  t.(e,L.)  *  -A  if  A  <  1,  and 


e+0 


dt . 


t  (e,L  . )  ®  -A  if  A  >  1.)  Then,  for  L  near  L.  and  i  >  1,  lim  — —  (e,L)  <  C. 

3+1  3  J  C-0 

Proof :  As  in  Lenina  4.2,  we  use  Lemnas  2.1  and  2.2  alternately  to  see  that  increasing  L 
Increases  -t^  for  each  i  >  1.  I 

We  return  to  the  proof  of  Theorem  4.2.  Lemma  4.3  implies  that  for  L  <  Lj  (resp. 

L  >  Lj),  t^tCjL)  >  t^(E,L^)  (resp.  t^EjL)  <  t^tEjL..))  for  all  e  sufficiently  small. 

Also,  by  hypothesis,  lim  y(-1,£,L)  «  -lim  t.(E,L)  (resp.  -lim  t.  (£,L))  for  A  <  1 
£-*■0  E-M3  J  E-'-O 

(resp.  A  >  1).  Hence,  for  fixed  &  small,  and  all  £  sufficiently  small, 

y(-1,£,L.  -  6)  <  A  <  y(-1,E,L^  +  6)  . 

To  finish  i)  it  suffices  to  show  that  for  |L  -  L^|  <  S,  3y/3L  >  0.  This  is  somewhat 
easier  than  the  analogous  proof  for  Theorem  4.1,  since  the  monotonicity  holds  for 
It  -  L  |  <  6,  where  o  may  be  chosen  independently  of  e.  (See  Figure  4.4.)  By 
definition  of  t^(£,L),  we  have 

y(ti(e,L),E,L)  *  -t^E,!)  . 

Differentiating  this  with  respect  to  L  yields 


3  dt- 

I*  (t^E.Ll.E.L)  =  -[1  +  y'(t  (e,t),E,L]  ~ 
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The  limiting  forms  of  solutions  to  (1.1)  with  the  same  B  but  different 
-L/e 

y'(0,£,L)  =  -e  (nearby  L). 

By  Lemma  4.2,  dt^/dL  <  0.  Also,  for  £,  6  sufficiently  small,  y* ( t  (e ,L) ,6 ,L)  is 

dti  1 

exponentially  small,  so  <3/3l) y(  t  ( e ,L)  ,e,L )  -  -  -jj—  <  0.  Furthermore,  as  e  *  0, 

y(t. (e,L) ,£,L)  +  lim  y(-1,e,L).  Hence  we  may  also  conclude  that  (3/3L)y(-1 ,e,L)  <  0. 

1  e-0  .  2  2 

ii)  To  show  that  (t^j  <  3B  «=>  (t^  )  <  3B  ,  it  suffices  to  show  that 

I tj  2 |  <  1 1^ I •  This  follows  from  the  arguments  in  the  proof  of  Theorem  4.1,  ii )  .  The 
assertion  for  A  <  1  follows  from  the  fact  that  the  first  point  to  the  left  of  t^  in  the 
solution  to  (3.1)  with  t^  =*  -A,  t^  =  -0  becomes  closer  to  t^  as  j  increases. 
Similarly,  if  A  *  1,  the  assertion  claimed  in  ii)  follows  from  the  fact  that  the  first 
point  to  the  right  of  t^+1  in  the  solution  to  (3.1)  with  tj+1  =  -A,  t,  >=  -B  is  closer 
to  tj+'j  aa  3  increases. 

iii)  This  follows  immediately  from  ii).  | 

The  two  combinations  left  to  be  done  (for  B  >  0  or,  if  B  <  0,  6  *>  0)  of  horizontal 
and  vertical  limiting  behavior  near  t  «  0,  t  »  -1  are  proved  by  the  methods  of  the  two 
previous  theorems.  Hence,  we  shall  merely  state  the  results.  Theorem  4.3  concerns 
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solutions  which,  in  the  limit,  are  horizontal  near  t  =  0  and  vertical  at  t  =  -1. 

Theorem  4.4  concerns  the  reverse  combination. 

t  1+1  > 

Theorem  4.1.  Suppose  3  >  0.  For  each  odd  j  >  3,  consider  the  solution  lt^  )  to 

(3.1)  with  tj+1  =  -1,  t^  =  -B.  Let  CV(B)  be  the  interval 

i)  There  is  a  unique  solution  to  (1.1),  (1.2)  having  a  boundary  layer  at  t  =  -1,  no 

boundary  layer  at  t  “  0,  and  j  strictly  interior  turning  points,  if  A  f  CV(B)  and 

2  2  2  2 
(tj)  >  33  .  A  necessary  condition  for  such  a  solution  is  A  e  cf  C  (B),  (t^)  >  3B  . 

ii)  If  A  *  1  or  B  >  0,  there  are  at  most  finitely  many  j  for  which  there  is 

such  a  solution.  I 

Theorem  4.4.  Suppose  B>0  or,  if  B  <  0,  6=0.  For  each  odd  j  >  3,  consider  the 
solution  {t?}  to  (3.1)  with  tj  =  -A,  tQ  =  0.  Let  ( A )  be  the  interval  (-“,-t^ ) . 

Then 

i)  There  is  a  solution  to  (1.1),  (1.2)  with  a  boundary  layer  at  t  ”  0,  no  boundary 

layer  at  t  =  -1  and  j  strictly  interior  crossings  with  y  =  -t  if  Be  A)  and 

tj+1  <  -1  (for  A  <  1)  or  -1  <  fcj  <  0  (for  A  >  1).  A  necessary  condition  for  such  a 

solution  is  B(  d  DAh ),  t^+1  <  -1  (for  A  <  1)  or  -1  <  tj  <  0  for  A  >  1. 

ii)  If  A  *  1  or  B  >  0,  there  are  at  most  finitely  many  j  for  which  there  is 

such  a  solution.  I 

The  techniques  of  Section  3  used  in  Theorems  4.3  and  4.4  work  only  for  j  >  3.  There 

are  solutions  for  j  ”  1,  and  these  cases  may  be  done  by  going  back  to  the  results  of 

Section  2  from  which  those  of  Section  3  were  derived: 

1  —  1  2 

Proposition  4.1.  Suppose  0  <  B  <  —  .  Let  L  «  —  -  B  ,  and  y  (t,6)  the  solution  to 

j-  13  1 

(1.1)  with  y  (0,c)  ”  B,  yj(0,e)  »  -1.  Let  £j(B)  be  the  open,  first  vertical  segment  of 

the  limiting  solutions  of  y(t,c)  as  e  *  0.  (See  Figure  4.5.)  Then  if  A  e  C^tB), 

there  is  a  unique  solution  to  (1.1),  (1.2)  with  a  boundary  layer  at  t  «  -1,  no  boundary 

layer  at  t  »  0,  and  intersecting  y  «  -t  exactly  once.  Furthermore,  if  y(t,e)  is  this 

-L1 (e)/c  _  _  1 

solution,  then  y'(0,c)  -  -e  ,  where  L,(c)  ♦  L,  as  e  +  0.  If  B  >  —  there  is 

1  1  /3 


no  such  solution 


Figure  4.5 


A  limiting  solution  for  the  case  j  =  1,  Theorem  4.3. 


Proposition  4.2.  Suppose  A  >  —  ,  B  <  A.  Then  there  is  a  unique  solution  to  (1.1),  (1.2) 

/3 

with  a  boundary  layer  at  t  *■  0,  no  boundary  layer  at  t  =  -1,  having  exactly  one  (resp. 

no)  strictly  interior  point  of  intersection  with  y  *  -t  if  A  <  1  (resp.  A  >  1),  and 

no  interior  points  of  intersection  with  y  =  t.  (See  Figure  4.6.)  Furthermore,  if  y(t,e) 

is  this  solution,  then  y'(0)  “  -K  !e)/e,  where  K  (e)  *  K  i  (A3  -  B3)/3.  If  A  <  — •  , 

/3 

there  is  no  such  solution.  | 


y 


A  limiting  solution  for  the  case  j  »  1,  Theorem  4.4. 


The  four  previous  theorems  Joel  with  solutions  to  (1.1)  for  which  there  are  no 

strictly  interior  crossing  points  with  y  =  t,  i.c.  B  >  0,  or  B  <  0  and  0  “  0.  We 

now  indicate  how  those  theorems  must  be  modified  if  B  <  0  and  6  <  0.  The  next  theorem 

refers  to  solutions  which,  in  the  limit  as  e  ♦  0 ,  have  a  vertical  segment  at  t  *=  0;  the 

final  one  deals  with  solutions  that  approach  a  horizontal  line  near  t  =  0. 

Theorem  4.5.  Suppose  B  <  0.  Consider  the  solution  (0,t^)  to  (3.1),  (3.9),  (3.10) 

with  tj+|  =  -1  it  the  limiting  behavior  near  t  =  -1  is  vertical,  or  t^  =  -A  (resp. 

tj+.j  »  -A)  if  the  limiting  behavior  near  t  =  -1  is  horizontal  and  A  <  1  (resp. 

A  >  1).  Let  i 4  be  the  interval  it  the  limiting  solution  is  vertical 

near  t  »  -1.  Let  3 »  (-*,6).  If  the  limiting  solution  is  horizontal  near  t  •  0,  let 
DAk)  *  (-*>,0).  Then 

i)  Theorems  4.1  i)  and  4.4  i)  continue  to  hold  if  A^t  B ^  and  Dj(A)  are  replaced 
by  A^,  B j  and  0_(A),  and  the  word  "odd"  changed  to  "even". 

ii)  V  J43  J6  ?2C?4C  —  • 

iii)  For  each  A  *  1,  there  are  finitely  many  j  for  which  A  e  ,3^.  | 

Theorem  4.6.  Suppose  B  <  0.  Consider  the  solution  {0,t?}  to  (3.1),  (3-9)  with  0  *  B 

and  tj+.j  or  t^  given  as  in  Theorem  4.5.  Let  Tt s  be  as  in  Theorem  4.5.  Then 
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i)  Theorems  4.2,  4.3  i)  continue  to  hold  if  A ^  is  replaced  by  Ay  (t2)2  >  3B2 
replaced  by  (t.j)2  >  3B2,  and  "odd"  is  replaced  by  "even". 

ii)  As  in  Theorem  4.5. 

iii)  As  in  Theorem  4.5.  I 

Remarks :  1.  For  each  of  the  four  combinations  of  horizontal  and  vertical  behavior  near 

the  endpoints  (6  cases  for  B  <  0)  there  is  a  slight  difference  between  necessary  and 
sufficient  conditions;  the  necessary  conditions  involve  certain  closed  intervals,  while  the 
sufficient  conditions  hold  on  the  interior  of  those  intervals.  The  difference  is  due  to 
the  fact  that  for  each  fixed  e  >  0,  there  is  a  region  around  certain  values  of  A  and 
B  for  which  the  number  of  solutions  is  not  necessarily  the  same  as  the  number  of  limiting 
solutions  as  c  *  0;  the  size  of  each  such  region  goes  to  0  with  c. 
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To  see  why  this  is  so,  it  is  useful  to  consider  (1.1)  as  an  initial  value  problem, 

with  B  =*  y(0,£)  and  y'(0,£)  specified.  As  an  example,  we  shall  consider  B  fixed, 

B  £  B.,  B  /  B,  (i.e.  B  <  - —  and  close  to  —  ),  and  increase  |y'(0,e)|  from  0. 

/3  _  /f 

For  y’  =  0,  y(t,c)  =  B.  Jf  y'(0,c)  «  -e  ,  and  L  is  now  decreased  from  •»,  there 
is  a  last  value  of  L  (namely  L  =  1/3  -  B2)  for  which  lim  y  =  B.  For  fixed  c  small, 

_  e*o 

y(-1»e,L)  stays  o(1)  from  y  *»  B  for  L  >  L. .  As  L  passes  ,  y(-1,c,L)  moves 
rapidly  up.  (By  "rapid",  we  mean  that  a  change  of  order  0(1)  in  y(-1,e,L)  is  produced 
by  a  o(1)  change  in  L.)  This  increase  continues  until  y(-1,£,L)  hits  some  point  close 
to  the  upper  endpoint  of  C.j(B),  and  y(-1,e,L>  starts  to  decrease.  Since  the  upper 
limit  is  not  exactly  the  boundary  of  C^(B)  for  fixed  c,  the  number  of  solutions  to 
(1.1),  (1.2)  for  A  very  close  to  this  boundary  may  be  two  more  or  less  than  the  number 
of  limiting  solutions. 

After  L  passes  ,  a  decrease  in  L  leads  to  a  slow  decrease  in  y(-1,£,L).  (By 

"slow",  we  mean  that  a  change  of  0(1)  in  y(-1,£,L)  is  produced  by  an  0(1)  change  in 

L.)  (See  Figure  4.7.)  The  decrease  continues  until  L  reaches  a  value  •  L  where  the 

-Z3/c 

solution  to  (1.1)  with  y(0,e,L)  =  B,  y'(0,£,L)  =  ~e  has  a  limiting  behavior  with  a 

vertical  segment  at  t  *  -1  and  three  strictly  interior  turning  points.  It  then  undergoes 
a  rapid  increase  in  y(-1,£,t)  for  L  «  1,^  and,  as  before,  stops  at  a  value  of  y  near 
the  upper  endpoint  of  Cg( B) .  This  alternation  of  rapid  increases  and  slow  declines  in 
y(-1,£,L)  continues  a  finite  number  of  times  for  L  e  [0,“].  (Indeed,  for  B  close 
enough  to  —  ,  the  first  rise  near  L  «  L  is  the  only  such  rise;  the  number  of  rises 

il  1 

can  be  computed  by  solving  (3.1)  with  tg  =  0,  t^  *=  B,  and  finding  the  last  j  such 
that  tjj  >  -1.)  The  rapid  rises  correspond  to  solutions  which  are  approximately  vertical 
near  t  «  -1;  the  slow  declines  correspond  to  solutions  which  are  close  to  horizontal 
near  t  •  - 1 . 
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4.7a  4.7b 

Figures  4.7a,  4.7b 

The  limiting  forms  of  solutions  with  y(0,£)  •  B,  y*  (0,e)  «  -e  L^e ,  with 
varying  L.  The  arrow  indicates  motion  as  L  is  decreased;  the  double 
arrow  indicates  rapid.  Figure  4.7a:  L  “  }  Figure  4.7b:  <  L  <  L ^ . 

Thereafter,  as  ]y'|  increases,  we  consider  y’(0,e)  =  -K/E.  Increasing  K  (see 
Figure  4.8a)  at  first  increases  y(-1,e,K)  slowly  (slow  and  rapid  have  the  same  meaning  as 
for  y(-1,e,L)).  The  slow  increase  ends  when  y(-J,e,K)  reaches  a  point  near  the  upper 
value  of  Aj,  at  K  ■  K  ;  for  K  near  K1#  there  is  then  a  rapid  decrease  in 
y(-1,e,K),  (see  Figure  4.8b)  which  ends  near  the  bottom  of  Aj.  For  larger  K,  the 
value  of  y(-1,e,K)  thereafter  increases  with  K,  and  is  unbounded.  If  B  t  4 ^ , 

B  f  ^j+2  *°r  some  j»  there  would  be  a  finite  sequence  of  slow  increases  and  rapid 

declines  in  y(-1,e,K)  before  the  final  monotone  increase.  Once  again,  near  every  value 
of  A  for  which  (3/3K)y(-1 ,c,X)  =■  0,  the  asymptotic  count  of  solutions  to  (1.1),  (1.2) 
nay  be  off  by  two. 

2.  Consider  (1.1),  (1.2)  now  as  a  boundary  value  problem,  with  A  and  B  as 
variables.  For  example,  suppose  we  fix  B  and  c  >  0  with  e  <<  1.  As  we  vary  A,  the 
number  of  solutions  changes  by  two  (either  up  or  down)  as  certain  critical  values  of  A 
are  passed;  these  points  correspond  to  changes  in  the  sign  of  (3/3K)y(-1,e,K)  or 
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(i/cL)y(-1 ,e,L)  and  have  been  computed  n  the  limit  as  e  *  0 .  The  disappearance  (or 
appearance)  of  two  solutions  is  accomplished  by  a  bifurcation:  two  solutions  coalesce  and 


then  disappear.  The  prototype  examnle  is  B  <  —  ,  A  variable.  As  A  decreases  past 
,  ‘  /3 

—  ,  the  left  boundary  layer  of  the  solution  which  is  asymptotically  i  A  (for 
/  3 

-1  <  t  <  0)  tends  to  zero  in  length,  and  the  solution  coalesces  with  another  solution 

which  is  asymptotically  =  —  (for  -1  <  t  <  0).  (See  Figure  4.9.)  For  A  <  •—  and  c 

/  3  /3 

small  enough,  neither  of  these  solutions  exists. 


Figures  4.8a,  4.8b 

The  limiting  forms  of  solutions  with  y(0,c)  «  B,  y* (0,e)  *  -K/e,  with 
varying  K.  The  arrows  indicate  motion  of  y(-1,e,K)  as  K  increases. 
Figure  4.8a:  K  <  K^;  Figure  4.8b:  K  »  • 


3.  The  methods  of  this  paper  can  also  be  used  to  get  a  complete  analysis  of  solutions 
for  t  >  0. 
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Appendix :  The  cape  B  <  A  <  0. 

Proposition  A.  Suppose  B  <  A  <  0.  Then  there  is  a  unique  solution  to  (1.1),  (1.2).  The 

limiting  behavior  cf  the  solution  for  -1  <  t  <  0  is  as  follows: 

if  B  <  -  <  A,  lim  y(t,£)  =  -  ~  # 

✓I  C-*0  /  3 

if  B  <  A  <  -  —  ,  lim  y(t,£)  =  A  , 

/3  e+0 

if  -  —  <  B  <  A,  lim  y(t,£)  =  B  . 

/3  e+0 


(See  Figure  A. 1 . ) 
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Figure  A. 1 

The  limiting  forms  of  solutions  for  the  three  different  cases  of 
Proposition  A. 

Proof:  It  follows  from  the  estimates  of  Section  2  that,  in  each  case,  the  above  constant 
is  the  only  possible  limiting  behavior.  To  show  that  there  is  a  unique  solution 
corresponding  to  these  boundary  conditions,  we  may  use  the  arguments  of  Theorem  4.1  and 
Theorem  4.2.  Alternatively,  we  may  apply  the  maximum  principle  to  the  difference  of 
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two  solutions.  If  yr  y2  are  solutions,  then  w  =  y1  -  y2  satisfies 

cw"  -  (y2  -  t2)w’  -  lyjtyj  +  y,))w  =  0  » 
w(-1  )  •>  w(0)  -  0  . 


Since  [y2<y2  +  y2>]  >  0  we  see  that 


w  =  0  • 
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